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ABSTRACT 

We present simultaneous UV-G-R-I monitoring of 19 M dwarfs that revealed a huge 
flare on the M9 dwarf with an amplitude in the UV of at least 6 magnitudes. This is 
one of the strongest detections ever of an optical flare on an M star and one of the first 
in an ultracool dwarf (UCD, spectral types later than about M7). Four intermediate 
strength flares (A?nuv < 4 mag) were found in this and three other targets. For the 
whole sample we deduce a flare probability of 0.013 (rate of 0.018hr -1 ), and 0.049 
(0.090 hr -1 ) for 2M1707+64 alone. Deviations of the flare emission from a blackbody is 
consistent with strong H n line emission. We also c onfirm the previously found rotation 
period for 2M1707+64 llRockenfeller et al ]m) and determine it more precisely to 
be 3.619 ±0.015 hr. 
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1 INTRODUCTION 

Numerous early and mid M dwarfs are known to show short 
timescale brightness eruptions flares similar to those on the 
sun. They are probably caused by magnetic fields, whereby 
energy is locally converted to thermal energy by a reconnec¬ 
tion of field lines, the consequence of which is a significant 
increase in the total stellar luminosity. Flares manifest them¬ 
selves in various wavelength regions, from the radio to the 
X-ray and as line or continuum emission. The corresponding 
motions of plas ma have been det ected by emission line asym¬ 
metry llFuhrmeister et a DiS e.g.). In the recent years, a 
few flares in the optical as well as emission line flares have 
also been observed in late M and L dwarfs (see below). 

The creation process of magnetic fields in fully convec¬ 
tive objects (such as UCDs) is not yet fully understood. But 
it must differ from the aR-dynamo present in higher mass 
stars because this depends on the interface between radia¬ 
tive and convective zones. However, the detection of flares as 
well as cool surface spots strongly suggests the existence of 
magnetic fields in UCDs. Recently, Chabrier & Kiiker (astro- 
ph/0510075) modelled an a 2 -dynamo for ultracool dwarfs. 
Among other things, they are able to explain the saturation 
in the rotatio n-activity relation which has been observed i n 
these objects iDelfosse et al.lll998i : iMohantv fc Basrill2003lf . 

Amplitudes of flares on K and M dwarfs in the optical 
have been observe d to increase strongly from the red (I - 
band) to the U/UV llStenanov et alll99.4lEason et alll992l) 


and no signatures can be traced in the near-infrared. The 
duration of these events seems to positively correlate with 
their amplitude and ranges from seconds to a few hours. 

Although several extensive photometric monitoring pro¬ 
grams of L and T dwarfs have been carried out, only 
a few flares have been reported, probably because most 
of these surveys were conducted in the I-band or near- 
infrared (JHK). Nevertheless, flares have recently been de¬ 
tected by optical broadband photometry, for example by 
iKoen (2005b) in an M8.5 d warf with an R-band amplitude 
of almo st 2 mag and by lScholz fc Eisloffell (j2005) in a young 
VLMS. E oe 3 ''20115a i reports of a possible flare event in an 
object as late as L8 (R-band amplitude of 0.1 5 mag). Emis¬ 
sion li ne flares have also been reported , e.g. b v ILiebert et alJ 
jl999|) (M9.5 object). iL iebert et, alJ feOO.'j) (L5) and by 
iFuhmieister fo Schmit,dT20fulW M9). 

iRockenfeller et aD ll200flf (henceforth R06) found peri¬ 
odic variability in the GRI-bands of 2M1707+64 at a pe¬ 
riod of 3.65 ± 0.1 hr with an amplitude of 0.014 mag in the 
I-band. By modelling the amplitudes they showed the mod¬ 
ulation was best explained by localised cool surface spots 
in this late-type object (which has a spectral type and es¬ 
timated age placing it near to the hy drogen burn ing limit). 
H a quiescent emission was found by iGizis et~aL feoQQ I) at 
an equivalent width of 9.8A, which suggests 2M1707+64 is 
a very active M dwarf, even during its non-flaring state. 

Here, we report the detection of two flares on 
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Figure 1. The total efficiencies of the BUSCA UV, G, R and I 
channels, taking into account the dichroics, the CCD efficiency 
and the Bessell-I filter used in the I-band. The efficiency of the 
UV filter drops off rapidly blueward of 0.3 /im due to the atmo¬ 
sphere. Over plotted for re f ere nee is a spectrum of the M4.5 dwarf 
DENIS P-J1158-1201 fromlMartin et alJ il99St) . 

2M1707+64 of which one has exceptionally large amplitudes 
throughout the optical, with a peak of more than 6 magni¬ 
tudes in the UV. We also report the detection of smaller 
(but still large) flares on three other late M dwarfs. 


2 DATA ACQUISITION, REDUCTION AND 

ANALYSIS 

Simultaneous multichannel UV,G,R,I observations of 19 
M dwarfs were obtained over 20 nights in June 2002 and 
2003 with the four channel CCD camera BLTSCA at the 2.2m 
telescope of the Calar Alto Observatory, Spain. BUSCA uses 
dichroics to split the light beam into the four wavelength 
bands which are imaged by four separate CCD cameras. 
Fig. □ shows the passbands. The UV, G and R passbands 
are defined by the CCD response and dichroic transmission 
function; for the I band we additionally used a Bessel I filter. 

Follow-up observations of 2M1707+64 were performed 
in June 2005 on two nights (14th and 15th) - again in the 
four channels with BUSCA - during Directors Discretionary 
Time, as well as on three nights (24th, 26th and 29th June 
2005) in the I-band with the 1.5 m telescope at Maidanak 
Observatory, Uzbekistan. All data were reduced in the same 
way as described in R06, except that the Maidanak data 
required an additional fringe subtraction. As in R06, light 
curves were obtained from differential photometry using a 
set of stable reference stars. Tabled shows the observation 
log. See R06 for detailed information on the analysis of the 
2002/03 BUSCA data, in particular time series analysis, de¬ 
termination of rotation periods and spot modelling. 


Figure 2. Relative light curves of 2M1707+64 on 2002.06.02, 
showing the I-band (solid line), R-band (dashed line) and G-band 
data (dotted line). The UV-band flare amplitude is larger than 6 
magnitudes. The error bars are also plotted, but except for the 
G-band are too small to be visible. 

3 FLARES ON 2M1707+64 

The time series of 2M1707+64 shows a large increase in 
brightness toward the end of observations on the morning 
of 2002.06.02 with increasing amplitudes from I to UV (see 
Fig- 13- The wavelength dependence, the shape of the light 
curves and the duration of the event (more than 1.0 hr) are 
consis tent with the signatures of stellar flares on late-type 
stars dStepanov et alll995T) . Since the time sampling is quite 
coarse with respect to the duration of the event, our obser¬ 
vations probably missed the maximum. Our reported ampli¬ 
tudes are therefore lower limits to the maximum flare am¬ 
plitudes. Furthermore, because the target is not visible in 
the UV-band during its quiescent state we have estimated a 
lower limit to the UV-amplitude by equating the quiescent 
target magnitude to the limiting magnitude of the exposure 
(which is approximately the magnitude of the faintest star 
visible). The very large amplitudes in all four channels with 
no change in the reference stars prove that the brightening is 
intrinsic to the source and is not a resul t of second ord er ex¬ 
tinct ion effects, which are far smaller dBailer-Jones fc Lamml 
120031) . The brightness evolution in the UV-band can be seen 
in Fig. [3 

The flare amplitudes in the I-, R-, G- and UV-band 
are 1.0, 2.7, 2.9 and 6 magnitudes, respectively, which 
is exceptionally large (compare the papers mentioned in 
the introduction). This is one of the strongest broad¬ 
band optical flares ever reported (more so as we probably 
missed th e maximum of th e eve nt). As a comparison, see 
iBatvrshinova fc Ibragimov d200lll who recorded a superflare 
on an M4 dwarf with an amplitude in the U-band of 11 mag. 

If we initially assume that the radiation from the flare 
is thermal and so obeys a blackbody distribution, we can 
attempt to determine the plasma temperature (T p ) and the 
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Table 1. Observation log of the four M flare dwarfs, listing the objects’ names, the dates of ob¬ 
servations, the duration of observation the individual nights (in hours) and the instrument/site. 


object 

dates 

duration [hr] 

instrument 

2MASSW J1707183+643933 

2002.06.01 / 02 

6.4, 5.6 

BUSCA 


2005.06.14 / 15 

5.7, 3.9 

BUSCA 


2005.06.24 / 26 / 29 

5.6, 5.1, 4.4 

Maidanak 

2MASSW J1344582+771551 

2002.05.29 / 2005.06.01 

5.2, 6.2 

BUSCA 

2MASSW J1546054+374946 

2003.06.06 / 07 /08 

6.4, 2.6, 2.4 

BUSCA 

2MASSW J1714523+301941 

2003.06.03 / 05 

6.1, 5.4 

BUSCA 
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Figure 3. Set of four UV-band images of 2M1707+64 on 
2002.06.02. The first image (upper left) shows the object during 
quiescent state, whereas the following three show the first three 
(out of four) images of the flare event — the numbers denote the 
time-order. Each image is a Ri 2.4' ■ 1.8' clip from the whole field, 
with North up and East to the left. 


(projected) fractional area of the stellar sur face that is ef¬ 
fected by the flare (A). We used equation 1 of lde .laser et all 
( 19861) for an optically thick plasma to simulate flare ampli¬ 
tudes over a wide parameter space (5000 K < T p < 10 5 K 
and 10- 7 < A < 0.5). Because of the low effective tem¬ 
perature of 2M1707+6 4 , we used an appropriate synthetic 
spectra of I Allard et al.. f2001 1 to model the flux during the 
quiescent state. The goodness-of-fit of an individual model 
was judged by a simple function f(T p ,A) = 1/4 • (|<5uv| + 
|<5g| + |<5r| + N), where 5 i is the difference in magnitudes be¬ 
tween the observed amplitude and the simulated amplitude 
in channel i. This reveals a degeneracy in the model (i.e. 
doubling the temperature and halving the area leads to the 
same value of /). Moreover, the large values of f{T p , A) for 
all fits indicate that the blackbody approximation for this 
flare’s emission is inappropriate (and/or that the quiescent 
spectral model is poor). Even the best-fitting simulations 
could only fit 2 out of 4 channels to within 0.1 mag (UV and 
G-band), the third within 0.6 mag (I-band) and the R-band 
always deviated by more than 1.8 mag. This could partly 
be explained by strong H a line emission, which falls within 
the R-band. However, one has also to consider that the UV- 
amplitude is poorly constrained (we only have a lower limit) 
which limits the validity of the fitting procedure. Depar¬ 
ture from bl ackbod y is no t unexpected. For example, from 
Table 2 of lde Jager et all (1986), we can similarly deduce 


non-blackbody emission for a flare on the K6V dwarf BY 
Draconis. 

The BUSCA 2005 data reveal another flare on 
2M1707+64 at the last data point on the morning of 
2005.06.14. For this second flare we find (lower limits to 
the) amplitudes of 3.70, 0.94, 0.77 and 0.16 magnitudes in 
the UV, G, R and I-bands respectively. Although at the 
very end of observations, the size of the amplitudes and the 
absence of variations in the reference stars rule out an ob¬ 
servational artifact. Again, the model fits with the lowest 
/-value show the R-band amplitude to deviate by more than 
0.5 mag while the other three bands can be reproduced to 
better than 0.1 mag. 

In total we obtained 22.3 hr of BUSCA data in 2002 and 
2005 on 2M1707+64 which revealed two flares. With a total 
recorded flare duration of 1.1 hr this implies a probability 
that the object is flaring at any one time to be 0.049. Note, 
however, that the flare continued after the end of observa¬ 
tions on that morning, so this is a lower limit (see Figl^Jl. 
If we instead just count flares, we deduce a flare rate of 
0.090 ± 0.063 hr -1 where the error is obtained by counting 
(Poisson) statistics. If we also take into account the 15.2 hr 
of Maidanak data whose sensitivity is restricted to strong 
flares (approx. Ami > 0.05, compare Table. we obtain a 
flare rate of 0.029 (0.053 ± 0.038 hr -1 ). 


4 PERIODIC VARIABILITY IN 2M1707+64 

From the five nights of follow-up observations in June 2005 
we detect periodic variations with periods (obtained from 
Scargle periodograms) in the individual nights ranging from 
3.61 to 3.64 hours at a 99.9 % confidence level. This con¬ 
firms the period we claimed in R06. The application of the 
CLEAN periodogram to the combined BUSCA-Maidanak 
data of 2005 yields a period of P = 3.619±0.020 (±0.015) hr. 
The first uncertainty estimate (0.020 hr) was obtained by 
Monte-Carlo simulations which account for the time sam¬ 
pling (see R06). The second estimate is from the theoret¬ 
ical expression for the width of peaks in a periodogram, 
A P = P 2 /(2T), where P is the period and T the total 
time span covered by the observations. This is only valid for 
uniform time sampling, an assumption which is strongly vi¬ 
olated here. The former result overcomes this, although still 
assumes the signal to be sinusoidal with Gaussian noise. 

We can exploit the non-uniform time sampling in this 
case to make a more robust estimate of the period uncer¬ 
tainty by constraining a single sinusoidal light curve to fit 
the data on all nights. The two consecutive BUSCA nights 
separated by several days from the three Maidanak nights in 
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2005 (see Tabled gives a strong constraint on this fit. The 
best fit gives a period of 3.619 hr, the same as found with 
the CLEAN periodogram. To estimate the uncertainty in 
this, we fix the phase of the sinewave at 2005.06.27 00:00UT 
(= to) and vary its period until we achieve a phase shift of 
7r over the 13 day difference between the 14/15 June ob¬ 
servations and to. The difference between this period and 
the nominal one is 0.015 hr; we (conservatively) define this 
to be the uncertainty. It cannot be much larger, because if 
we further change the period until we reach a phase shift 
of 27 t (which still gives a reasonable fit on 2005.06.14 & 
15), we now see a significant deviation on 2005.06.24 & 
29. Hence, under the assumption that the observed peri¬ 
odic variability on all five nights in June 2005 was caused 
by the sa me star spot(s), a nd that differential rotation is 
negligible llBarnes et aljEooiill . we establish the rotation pe¬ 
riod (of the spot system) to be 3.619 ±0.015 hr, an accuracy 
of better than 1 min. Although the periodogram implies a 
similar uncertainty, the sinewave fit is more robust. 1 The 
lifetim e of ste llar surface sp o ts ca n be much longer than two 
weeks llStrassmeier fc Bopdll99^l so our assumption is rea¬ 
sonable. On the other hand, the lifetime of spots on very low 
mass stars and brown dwarfs has been little studied and may 
differ from G and K stars. If the spot distribution evolves 
significantly during the 2005 observations then the fitting 
method is not valid. The data from the five nights phased 
to 3.619 hr are shown in Fig. 2] 
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Figure 4. 2M1707+64: I-band light curve of the five nights be¬ 
tween 2005.06.14 and 29 phased to a period of 3.619 hr. A typical 
error bar is also shown. Note that two cycles of the phased light 
curve are plotted. 


5 FLARES ON OTHER OBJECTS 

From our 20 night BUSCA survey (R06), we found three fur¬ 
ther flares, one in each of the following objects: 2M1344±77, 
2M1546±37 and 2M1714±30. All three dwarfs have de¬ 
tectable quiescent H a emission dGizis et al . 20 q 3). Ampli¬ 
tudes and further properties are listed in Table [5] R06 
detected variability in the R-band in 2M1546±37 and 
2M1714±30 at an amplitude of 0.012 mag in both cases. A 
tentative period of 6.9 hr was assigned to the latter in the R- 
band only. In all three cases the blackbody model fits the ob¬ 
served amplitudes relatively well and the following applies to 
all three flare events. If we again consider those simulations 
with low /-values, we find two classes of results. One well 
reproduces I, R and G amplitudes with a strong deviation 
in the UV-band (« 1.1 mag) and the other fits the I, R and 
UV well with a deviating G-band amplitude (« 0.2mag). 
Since the UV-amplitude is not as well determined as the 
other three (see Sect. CD, the former seems to be more likely 
and thus indicates relatively good agreement with a black- 
body emission. However, various emission lines have been 
observed in flare spectra in both bands llEason et al ] iH) 
and could thus modify the continuum emission (that itself 
also could be non-blackbody) such that it deviates from a 
blackbody. This could explain any residual deviations. 

1 Mathematically the uncertainty derived from the sinewave¬ 
fitting is “ X P, identical to the theoretical expression for the 
periodogram. However, the derivation of the periodogram expres¬ 
sion depends on the (invalid) assumption of uniform sampling. 
With the fitting method, in contrast, we explicitly determine the 
quality of fit of the data to the model at individual epochs for 
arbitrary phase shifts. 


6 CONCLUSIONS 

From 227.6 hours of multiband monitoring of 19 M dwarfs 
(218.0 hr with BUSCA in 2002 (R06) and 9.6 hr with 
BUSCA in 2005), we detected five flares on 4 different tar¬ 
gets. From this we infer a flare probability at any one time of 
0.013 and a flare rate of 0.022 ± 0.010 hr -1 for field M2-M9 
dwarfs. Although the observation time per object is limited 
(less than 14 hours per object, except for the follow up ob¬ 
servations of 2M1707±64), we conclude that flare activity is 
not rare in old (field) mid to late M dwarfs and might even 
be present in most of these objects. 

We discovered a massive flare on the M9 dwarf 
2M1707±64 with a UV amplitude of more than 6 mag. This 
object, which is either a very low mass star or brown dwarf, 
shows a repeat of flare activity at a rate of 0.090±0.063 hr -1 
(probability of flaring 0.049). With both this and another 
flare on the same object we found deviations of the flares’ 
4-band emission from that expected for a blackbody. In par¬ 
ticular, the R-band always fit poorly, which is consistent 
with the presence of strong H a line emission. 

We also improved the determination of the rotation pe¬ 
riod of 2M1707±64, deriving a value of 3.619 ±0.015 hr, un¬ 
der the assumption that the observed spot was stable over 
16 days of observations. 

We found three further flares on the late M dwarfs 
2M1344±77, 2M1546±37 and 2M1714±30. Their ampli¬ 
tudes and energy release were much smaller than in the case 
of 2M1707±64 and their emission closer resemble that of a 
blackbody. The occurrence of flares and cool surface spots 
strongly suggests the existence of magnetic fields in UCDs, 
even though the exact creation mechanism remains enig¬ 
matic. 
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Table 2. Properties of the detected flares: spectral type, quiescent H a equivalent widths, the flares’ amplitudes in 
the four bands, an estim ate to their duration and the number of recorded data points (samples) for each flare. H a 
EWs are taken from lGizis et all |(2000j). 


object 

year 

SpT 

H a EW [A] 

A UV 

Aq 

Ar 

Ai 

duration [hr] 

samples 

2M1707+64 

2002 

M9 

9.8 

> 6 

2.9 

2.7 

1.0 

> 1.0 

4 

2M1707+64 

2005 



> 3.7 

0.94 

0.77 

0.16 

> 0.1 

1 

2M1344+77 

2002 

M7 

2.7 

> 1.6 

0.094 

0.047 

0.007 

< 0.5 

1 

2M1546+37 

2003 

M7.5 

10.9 

> 1.7 

0.155 

0.121 

0.022 

< 0.75 

1 

2M1714+30 

2003 

M6.5 

5.4 

> 1.4 

0.075 

0.031 

0.007 

< 0.5 

1 
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